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Abstract 

A relativistic quantum mechanics is studied for bound hadronic sys- 
tems in the framework of the Point Form Relativistic Hamiltonian 
Dynamics. Negative energy states are introduced taking into account 
the restrictions imposed by a correct definition of the Poincare group 
generators. We obtain nonpathological, manifestly covariant wave 
equations that dynamically contain the contributions of the negative 
energy states. Auxiliary negative energy states are also introduced, 
specially for studying the interactions of the hadronic systems with 
external probes. 



PACS number(s): ll.30.Cp, 24.10.Jv, 03.65.Pm 



1. Introduction 

The study of hadronic few-body systems in terms of constituent particles rep- 
resents a very important tool for the determination of their intrinsic proper- 
ties and their interactions with external probes. We use the term constituent 
particle to mean a system that transforms as an irreducible representation of 
the Poincare group with definite mass, spin and internal symmetry quantum 
number. 

In this respect we recall that historically the investigation started with the 
study of light nuclei in terms of constituent nucleons. Later, many efforts 
have been devoted to the study of the hadrons in terms of constituent quarks. 
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For both cases nonrelativistic models were initially considered, later on, rel- 
ativistic corrections were added, improving the reproduction of the experi- 
mental data. 

Nowadays, the construction of hadronic covariant constituent models must 
be considered strictly necessary for an accurate description of these systems 
and for the study of their interactions with electroweak probes. In partic- 
ular, the present work, that represents a generalization of the Relativistic 
Hamiltonian Dynamics, is focussed on the study of quark (and quark-di- 
quark) models, but its formal developments can be also applied to the study 
of few-body nuclear systems. 

As for the relativistic covariant quark models, they represent effective models 
that should be related to quantum chromodynamics (QCD), incorporating 
in a nonperturbative way its symmetries and dynamical features. On the 
other hand, the nuclear models rely on a phenomenological meson exchange 
theory. 

1.1 General remarks 

Prom a theoretical point of wiew, we note that, for building few-body rela- 
tivistic models, two slightly different approaches can be followed. 

The first one, is denominated Relativistic Hamiltonian Dynamics (RHD). 
The second one is represented by the Relativistic Wave Equations (RWE). 

The aim of the first approach, that is the RHD, is to satisfy the Poincare 
Group commutation rules by defining the generators of that group in terms 

of the constituent particle operators. 

We now syntetically recall some technical aspects related to RHD. In the 
case of free particles the total generators of the Poincare group are given 
by the sum of the single particle generators. The main problem of RHD, 
not found in the nonrelativistic case, consists in fulfilling the Poincare group 
commTitation rules when the interaction is introduced in the generators. 
In this respect, RHD can be formulated in different ways, three of which have 
been introduced in the pionieristic work by Dirac [1]. They are called the 
Instant Form (IF), the Front Form (FF) and the Point Form (PF). 
The most relevant difference among them is represented by the way in which 
the generators depend on the interaction. 

In the IF both the Hamiltonian, i.e. the time translation generator, and 
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the boost generators are modified (with respect to the free case) by the 
interaction [2-4]. 

In the FF, hnear combinations of the components of the four-vectors are 
considered and, in consequence, not standard Lorentz transformations must 
be introduced. In this form of RHD, the interaction modifies pieces of both 
the boost and of the four-momentum [3] . 

Finally, in the PF [5-7], that is the scheeme adopted in the present work, 
the interaction modifies the total four-momentum, i.e. the Hamiltonian and 
the three-momentum of the system, that are given not only by the sum of 
the four-momenta of the costituents, but also receive a contribution that 
depends on the interaction operator. On the other hand, the boost is left 
free from interaction. Due to this last property, PF RHD has been defined as 
manifestly covariant [5]. This point will be analyzed in the next subsect.1.2, 
considering the definition of manifest covariance adopted in the present work. 

Theoretically, the relevant aspect of RHD is that its quantum mechanical 

properties are well defined, in analogy to nonrelativistic case. 
All the generators of the Poincare algebra are represented by hermitic op- 
erators so that the corresponding transformations are performed by unitary 
operators, satisfying Poincare covariance and allowing, at the same time, to 
use the standard quantum mechanical procedures for the calculation of phys- 
ical observables. The interaction operator is given by a quasipotential that 
is, in general, momentum dependent. 

Note that some aspects that are typical of quantum field theories have been 
completely excluded in RHD: in particular, the possibility of creating or de- 
stroying particles and the presence of negative energy states in the interaction 
amplitudes. In this work we shall analyze and overcome this second problem. 

The quark models based on RHD reproduce some general hadronic features 
related to QCD, like colour global symmetry, confinement and isospin invari- 
ance for the u d quark sector. 

Furthermore, in all the three forms of RHD very encouraging quantitative 
calculations have been done and are still in progress both for the hadronic 
spectra and for the electroweak form factors [8] . 

On the other hand, the dynamics of these models is intrinsically phenomeno- 
logical and partially unsatisfactory, because, as discussed above, some funda- 
mental aspects of QCD considered as field theory, are not taken into account 
by standard RHD. Moreover, the parameters that appear in the mass oper- 
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ator of these models are usually fitted to the experimental data in order to 

reproduce the hadronic spectra. 

The electroweak form factors are calculated in the relativistic impulse (or 
spectator) approximation by using a current operator that, in the case of 
the electromagnetic interaction, is not dynamically conserved. The discrep- 
ancies that remain between the theoretical calculations and the form factor 
experimental data are usually cured by inserting phenomenological quark 
form factors and/or vector meson exchange factors at the quark interaction 
vertex. 

The second possible approach for the relativistic study of few-body bound 
hadronic systems, is represented by the use of integro-differencial three- 
dimensional RWE. 

We leave aside from our discussion the four- dimensional Bethe-Salpeter equa- 
tion because, if it is not reduced to a three-dimensional form, its formalism 
cannot be easily interpreted in a physical way and also the practical solution 

of the equation presents serious difficulties. 

In general, the procedure to obtain the RWE (for two-body interacting sys- 
tems) starts from an underlying field theory. Historically, quantum electro- 
dynamics was considered. 

Prom the perturbative expansion of the field theory, a set of relevant Feynman 

graphs is selected, excluding, in general, the contributions due to the poles of 
the bosonic propagators. Usually, the box and (in some cases) the crossed box 
graphs are taken. The infinite series of these graphs, denoted respectively as 
ladder and crossed ladder, is summed up by means of a Lippman-Schwinger 
equation for the scattering matrix, in which a suitable Green function is used. 
Finally, from that equation, an integro-differencial, three-dimensional RWE 
for the bound states is derived. 

The origin of the differences among the various RWE is due to which Feyn- 
man graphs are selected and to the approximations done to sum them up 
in a three-dimensional form. The procedure outlined above establishes the 
connection between the RWE relativistic model and the underlying field the- 
ory. Generally, the solutions of the RWE cannot be interpreted as standard 
quantum mechanical wave functions. Also, discussing some specific examples 
we shall highhght some difficulties of the RWE approach. 

Historically, RWE have been introduced to study electromagnetically bound 
systems, like atoms and positronium, taking quantum electrodynamics as the 
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fundamental theory. Later, they have been applyed to the study of few-body 
nuclear systems. In this case it was used a phenomenological interaction 
operator that represents the exchange of pions and heavier mesons. 

The simplest case of RWE is represented by the positive energy states Black- 
enbeclcr-Sugar equation [9]. In the remainder of this section we shall not 
further discuss this case that is not specifically relevant for the main objective 
of this work. The Blackenbecler-Sugar equation will be considered again in 
sect. 4 (as a special example with no negative energy states) showing that it 
can be directly interpreted by means of the standard PF RHD formalism and 
re-written in a manifestly covariant way. 

Considering the Dirac-likc equation [10], we recall that, in this theo- 
retical difficulty was found because this equation presents an unphysical pole 
for a vanishing value of the total mass of the system. As a consequence, 
a correct normalization of the wave function is not possible. This problem 
has been analyzed and denoted as continuum dissolution and cockroach nest 
in refs. [11] and [12], respectively. An analysis of this problem will be also 
performed in the present work by according to the requirements of PF RHD. 
The Gross equation [13] represents a relevant example of a manifestly co- 
variant approach for the study of interacting systems. It has been obtained 
by means of a covariant procedure starting from the box and crossed box 
Feynman graphs. Even though also in this equation a pole for a vanishing 
total mass is present, this singularity is easily removed requiring the (phe- 
nomenological) kernel to vanish faster than the total mass. In any case the 
solutions of the equation are normalizable if the total mass is not equal to 
zero. The Gross equation is not manifestly symmetric, as such, under particle 
interchange, but it has been symmetrized in all its applications. Within this 
approach, the mass of the bound system is obtained as a pole below threshold 
of the interacting scattering matrix. 

We also consider the so-called Breit equation in the form originally given 
for the study of electromagnetically bound systems [14]. This equation is 
obtained from the four-dimensional Bethe-Salpcter equation approximating 
the electromagnetic interaction with an instantaneous quasipotential. Only 
box Feynman graphs are taken. This equation can be directly interpreted in 
the scheeme of RHD that will be developed in the present work. 
Finally, the Mandelzweig- Wallace integro-differential equation [11] is struc- 
turally similar to the previous one, but also includes, in the eikonal ap- 
proximation, the crossed box graphs, improving the correspondence with the 
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underlying field theory and obtaining the so called one-body limit when the 
mass of one particle is set equal to infinity. This equation can be also written 

in a covariant form [15, 16]. 

However, a standard definition of the Poincare group generators in terms of 
PF RHD quantum mechanical operators is not directly feasible within this 
model. Such problem arises because the Mandelzweig- Wallace equation is 
not an eigenvalue equation for the mass of the system. In other words, the 
mass operator is not defined explicitly. 

More recently, the Mandelzweig- Wallace formalism has been revised and de- 
noted as equal time reduction. In a clever work [17] concerning electron scat- 
tering on the deuteron, considered as a two-nucleon relativistic bound system, 
that equation has been written in a Hamiltonian form and the boost gener- 
ator is constructed, with some approximations, in the instant form RHD. 
With respect to this problem, in sect. 4 we shall show that, with some han- 
dling, an equation equivalent to the Mandelzweig- Wallace one can be written 
in the framework of our PF RHD model, otaining exact relativistic covari- 
ance. By means of this procedure we shall introduce the auxiliary states that 
represent another original aspect of our work. 

1.2 Motivations and structure of the work 

Having discussed the theoretical framework of PF RHD and RWE, we intro- 
duce the motivations of the present study. 

The long term objective of the relativistic quark model investigation would 
be the construction of a covariant Hamiltonian model for the study of the 
hadronic systems. This model should reproduce, with a limited number of 
parameters, the hadronic spectra and the response of the hadronic particles 
to clectroweak probes, that is elastic and inelastic form factors, Compton 
scattering amplitude, etc.. 

Relativistic covariance, possibly in manifest form, should represent the formal 
framework of the model. 

As for the dynamics, we point out that QCD is assumed to be the ultimate 
physical description for these system. For this reason, the covariant Hamil- 
tonian model should represent a solvable approximation of QCD in the sense 
that its (generalized) wave equation should be able to sum up the relevant 
QCD graphs, allowing to treat perturbatively other effects not included in 
the sum. 

The connection with QCD should allow to relate the parameters of the model 
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with those of QCD. In particular, we refer to the effective quark masses. Also, 
the quark- (anti) quark interaction should be possibly derived from QCD. 
These objectives are, obviously, ambitious and lie beyond the present under- 
standing of nonperturbative aspects of field theories. Deep and long investi- 
gations are required. However, the present work, as a preliminary step, can 
help to solve some specific theoretical problems, also improving the under- 
standing of the hadronic phenomenology. 

As for the main objective of the present work, that is focussed on the study of 
spin 1/2, identical, interacting quarks, we want to include, with some approx- 
imations, the negative energy states in PF RHD, in order to represent the 
dynamical rclativistic effects that have not been taken into account by stan- 
dard RHD models. We recall that in quantum field theories, like QCD, when 
the perturbative series is considered, the negative energy states appear in the 
propagators of the intermediate particles. For this reason, the procedure of 
the present work can help to improve the understanding of the link between 
RHD and the underlying field theory. Also in the case of few-body nuclear 
systems, the discussion about RWE shows that the contributions of the neg- 
ative energy states are considered necessary for constructing a dynamically 
consistent model. 

Technically, the insertion of negative energy states in the mass operator of the 
model, is obtained, in sects. 2 and 3, by generalizing the PF RHD construc- 
tion of the Poincare group generators, fulfilling, also in presence of negative 
energy states, the commutation rules of the algebra. As shown in sect. 3, 
this procedure is possible if all the particles of the state have the same en- 
ergy sign. These states represent the dynamical states of our generalized PF 
RHD. Such condition is strictly necessary in order to define the four- velocity 
operator V"^, that, in turn, is essential to introduce the total four-momentum 
of the system P>^. 

Our wave equation, written in sect. 4, is a coupled equation that involves 
positive and negative energy dynamical states and explicitly satisfies charge 
conjugation invariance. 

Another original result of this work is the introduction of the auxiliary states 
that represent the states in which particles with different energy signs are 
present. For this reason (as discussed before) they do not directly participate 
in the dynamics of the model. The auxiliary states are covariantly defined. 
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by means of the interaction operator, in eqs.(4.14b) and (4.18b) for two and 

three-body systems, respectively. 

In this way, the definition of the auxiliary states does not require the in- 
troduction of new parameters. Their expression is derived by using RWE 
formahsm as a link with field theory. 

In order to understand in more detail the physical meaning of the auxil- 
iary states, we recall that, in perturbative expansions, the negative energy 
terms of the fermionic propagators give rise to the so-called z-graphs, corre- 
sponding to intermediate states with one quark and a quark- antiquark pair 
[17-19]. The contributions of such states to the electroweak currents have 
been extensively studied in the nonrelativistic limit, within a constituent 
chiral quark model. In particular, the two-body (or exchange) electromag- 
netic four-currents have been derived [18], obtaining significant contributions 
for the inelastic transition amplitudes of all nucleon resonances. 
In the same context it has been also shown that the exchange contributions 
are necessary to satisfy Partial Conservation of Axial Current (PCAC) con- 
dition for the weak interactions, also leading to a possible interpretation of 
the missing nucleon spin as angular momentum carried by the nonvalence 
degrees of freedom of the nucleon [19]. 

Note that, while the electromagnetic exchange currents could be in part also 
derived by means of a minimal coupling substitution in the nonrelativistic 
Hamiltonian, on the other hand the derivation of the weak (axial) exchange 
currents strictly requires the use of the z-graphs. 

For the reasons discussed above, in our generahzed PF RHD the auxiliary 
states, that represent the pair terms in effective way, are expected to give rel- 
evant contributions to the interactions of the hadronic system with external 
probes. In subsequent works it will be studied at numerical level their effect 
on electroweak form factors. 

Furthermore, as for the possible relevance of both dynamical and auxiliary 
negative energy states for the electromagnetic interactions of hadronic sys- 
tems, we point out that Compton scattering amplitude, if studied by means 
of the Feynmam graphs of quantum electrodynamics, requires the presence 
of the negative energy states in the fermionic propagators, in order to obtain 
the correct expansion up to the second order in photon energy [20] . 
For this reason, photon scattering on nucleons (considered as bound quark 
states) can represent a very helpful tool for the study of some nontrivial 
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aspects of the hadronic dynamics that should be reproduced by RHD quark 
models. 

As for the strong interaction of the quarks, the negative energy states appear 
in the covariant wave equation of the model, given in eq.(4.6), at two levels. 
First, as discussed previously, we have the dynamical negative energy states. 
Furthermore, the connection established with RWE and, in turn, with the 
underlying field theory, gives rise to quadratic terms in the effective inter- 
action operators of eqs.(4.13b) and (4.17b) for two and three-body systems, 
respectively. These terms reproduce, in the wave equation, the effects due 
to intermediate states (of the scattering matrix expansion) with different en- 
ergy signs. By means of the same arguments used for the auxiliary states, 
the quadratic terms, representing the z-graphs [17], can be interpreted as 
contributions of quark- antiquark pairs. 

Note that (as in the case of the auxilary states) these (extra) quadratic terms 
do not introduce new parameters, being related to the linear terms of the 
quasipotential. The effects of the quadratic terms on the constituent inter- 
action have been studied in the context of few-body nuclear physics [13, 21]. 
In the case of quark models, the small value of the quark masses can give 
rise to highly nonrelativistic effects that were not present in the interaction 
of the nuclear systems. 

The kinetic term of our wave equation given in eq.(4.6) is similar to that of 
the Breit equation of ref.[14], where a perturbative technique has been also 
studied for the case of positronium, that is bound by the electromagnetic 
interaction. Taking into account that the properties of the wave equation 
are strictly related to the structure of the interaction operator, we note that 
in the present case, that is quark models, such operator is not directly known 
and should be determined taking into account the underlying field theory. 
The total effective interaction (linear -|- quadratic terms) is in any case highly 
momentum dependent and a carefuU formal and numerical examination of the 
wave equation should be performed to highlight its properties in connection 
with nonperturbative hadronic phenomenology. 

Another objective of this work is to write the wave equations of our gener- 
alized PF RHD in a manifestly covariant way. We define here as manifest 
covariance the property of an equation of being written in terms of quantities 
that (a) transform as tensors under Lorentz transformations and (b) are not 
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related to a specific reference frame. 

In this sense, tfie Bethe-Salpeter and Gross equations are manifestly co- 
variant. On the other hand, the wave equations of PF RHD, written in the 
standard formalism, do not fully satisfy the requirement (b) due to the use 
of the zero-momentum (rest) frame for the definition of the velocity states 
[5]. 

However, it is well-known that, if a theory is really covariant, with some 
handling it can be written in a manifestly covariant form, as it will be done 
for PF RHD in the present work. 

Obviously, manifest covariance only represents a formal property of the 
equations and does not lead, as such, to improve the knowledge of the 
hadronic dynamics. 

Technically, manifest covariance is obtained by means of (a) the definition of 
the projection states of eq.(3.9), where the three-momenta of — 1 particles 
and spatial part of the four-velocity of the system are selected as spatial 
variables, and (b) the choice of the normalization of the wave function given 
in eq.(4.5), leading as result to a great clarification the formalism of PF RHD. 
Note that the structure of the covariant integration in eq.(4.6) is similar to 
that originally introduced by the Gross equation [13]. 

We recall that the methods of PF RHD allow to define a consistent framework 
for the study of the electroweak interactions of the hadronic systems (elastic 
and inelastic form factors), specially by introducing as a starting point the so 
called relativistic impulse (or spectator) approximation [5]. Our manifestly 
covariant formalism also allows to calculate in a much more direct and clear 
way the matrix elements of the electroweak current in that approximation, as 
done in ref. [22] for a model with only positive energy states. Furthermore, in 
this way it is possible to introduce a minimal coupling procedure to derive a 
conserved electromagnetic current [23] for a model containing positive energy 
states. 

The formal developments of the paper are organized as follows. 

In sect. 2, we define the Poincare Group generators for the case of free particles 
also introducing the negative energy states. This objective is reached consid- 
ering that in the three- momentum Lorentz transformation (and. in general, 
in the Lorentz boost operators), the relevant parameter is the product Av, 
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introduced in eq.(2.2b), being A the energy sign of the particle and v the 
spatial part of the four-velocity boost parameter. 

We also study the connection between the state representation of PF RHD 
with that given by the standard Dirac equation spinors, which is currently 
used in the developments of quantum field theory. In particular we show 
the equivalence of the use of the RHD Wigner rotations with standard Dirac 
boosts for the calculation of the relevant matrix elements of the model. 

In sect. 3, we construct the generators of the Poincare algebra for interacting 
particles with negative energy states, verifying that the commutation rules of 
the Poincare Group are still fulfilled. To this aim, we introduce the definition 
of dynamical states and auxiliary states. Only the former enter in the wave 
equation of the model. We also introduce the projection states of eq.(3.9). 
These states depend on the spatial variables that are used to obtain the man- 
ifestly covariant form of the model. Furthermore, the interaction operator 
is written in terms of Dirac spinors and matrices highlighting its covariant 
character and allowing to derive its expression from field theories. 

In subsect.4.1, we write in eq.(4.6) the wave equation of the model in a man- 
ifestly covariant way. We show that relativistic covariance, realized by the 
Poincare Group commutation rules, as a dynamical consequence, automati- 
cally avoids, in our formalism, the continuum dissolution disease. 
In subsect.4.2 we use the techniques of the RWE (referring to the Mandelzweig- 
Wallace model) to establish a link with field theories. In particular we intro- 
duce, starting from the same interaction operator, both the auxiliary states 
in eqs.(4.14b) and (4.18b), and the quadratic terms in eqs.(4.13b) and (4.17b) 
for the effective quasipotential of the wave equation. 

In the appendix we also show explicitly that the interaction operator is in- 
variant under Lorentz transformations. 

2. Poincare group transfomations for single particle states 

As a starting point we recall the form of a Lorentz transformation (boost) 
of a four- vector — (/°, f). In all this work we consider canonical boosts. 
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The boosted four- vector has the form 



/° = ^;T + vf (2.1a) 

fft = f + v(vf ^ + /°) (2.1fc) 

These equations can be summarized as follows 

f,^ = L^Mr (2.1c) 

In the previous equations we have introduced the relative four-velocity — 
(v°, v) that satisfies VnV^ — 1 and, in consequence 

^;0 = ^;0(v) = [l + v2]V' i^-ld) 

The well-known relations v = 7U and — ^ — [1 — u^]"-*^/^ easily connect 
the four-velocity to the standard physical velocity u of the initial reference 
frame measured from the boosted one. 

We now consider the Lorentz transformation for the on-shell four-momentum 
of a single particle. In the quantum-mechanical model that will be studied in 
the following, the particle three-momentum p will be used as spatial variable. 
In order to study also negative energy states, we introduce the notation 
p^{X, p) = (Ae(p), p) being A the energy sign of the state and 

6(p) = [p^ + mT/^ 

the absolute value of the energy; also, m represents the mass of the particle. 
We enphasize that in our model the free energy of a particle p°(A, p) is always 
considered as a function of the energy sign A and of the three-momentum p 
of the state. According to eqs.(2.1 a,b) the Lorentz transformation is 

p°(A, p; v) = ^°Ae(p) + vp = Ae(p6(p; Av)) (2.2a) 

Pb(A, p; v) = pb(p; Av) = p + Av(Avp^^ + e(p)) (2.26) 
where v^, that is given by eq.(2.1d), takes the same value for A = +1, —1. 

Being the particle on-shell, the independent transformation is that of the 
three-momentum p, displayed in eq.(2.2b). For the development of the model 
we highlight the two following properties of this transformation : 



12 



(i) the transformation is not linear with respect to this variable p due to the 
presence of e(p) in the r.h.s of that equation; 

(ii) the transfomation only depends on the product Av. 

For the construction of the PF RHD we shall introduce the infinitesimal 
generator of the Lorentz boost, denoted in the present work as t. To this aim 
one has, first, to expand eq.(2.2b) up to the first order in the physical velocity 
u and, second, to construct the operator t that transforms the momentum 
eigenstates according to that expansion [1-4] . 

Furthermore, we point out that we shall construct a unitary representation 
of the boost generator in order to represent the Lorentz transformations 
according to the standard rule of quantum mechanics, i.e. analogously to 
nonrelativistic quantum mechanics. 

This procedure seems suitable for the study of systems composed of particles 
whose (strong) interaction is described by means of a quasipotential operator. 
On the other hand, we recall that the spinors of the Dirac equation are 
transformed by pseudounitary operators, that will be shown explicitly in 
eq.(2.22). In the Dirac equation formalism, unitarity is recovered in the 
context of field theory [24], whose construction is beyond the scope of the 
present work. 

However, we shall show that the two representations are connected by the 
transformation of eq.(2.18). Finally, we remark that the price paid for using a 
unitary representation of the boost operator is the nonlocal character of both 
the Hamiltonian operator and of the boost generator, that will be explicitly 
defined in eqs. (2.4 b-e) for positive energy states and, in general, in eqs.(2.10) 
and (2.11). 

We now introduce for a single particle all the 10 infinitesimal generators of the 
Poincare group, collectively denoted as [g^] (1=1,. ..,10). In more detail, these 

generators are: the 3-momentum p, i.e. the generator of the spatial transla- 
tion; the angular momentum j , i.e. the generator of the spatial rotations; t, 
i.e. the generator of the Lorentz transformations; finally, the Hamiltonian h, 
i.e. the generator of the time translations. Their commutation rules are 

[p-,p^] = [p-,h] = [r,h]=0 (2.3a) 
[J^/] = ^6"'^V (2.36) 
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=ie°/37p7 (2.3c) 

[r,p'^]=i5°''^h (2.3e) 

= -ie°/37^-7 (2.3/) 

By introducing the operator r canonically conjugated to p, that is 
[p^jT'^] = —i6°'^, and the spin operator s, the generators j, t and h can 
be put, for the positive energy states (A = +1), in the following form that 
satisfies the commutation rules of eqs.(2.3a-g). 

j = r X p + s (2.4a) 

h = e(p) (2.46) 

t = d + g (2.4c) 

with 

d=i(re(p)+e(p)r) {2Ad) 

- P ' (2.4e) 



e(p) + m 

As anticipated, the boost is represented by a unitary operator that, for a 
finite transformation, has the following form 



B{v) = exp(itw(v)) ~ 1 + itu (2.5) 

with 

w(v) = ^ tanh-^(4[) = A tanh-Vlul) (2.6) 
|v| \u\ 

Considering a state of three-momentum p, positive energy (A = +1) and 
z-projection of the spin a, the action of the boost operator on such state is 



S(v) |p, A = +1, <7 >= P^^^^V/'R{p; v) |p,(p; v), A = +1, a > (2.7) 

For these states we adopt the delta normalization shown below in eq.(2.8). 
In consequence, the numerical factor in the r.h.s. of the previous equation. 
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whose origin is due to the nonhnearity of the Lorentz tranformation with re- 
spect to p (see eq.(2.2a)), guarantees the correct normahzation of the boosted 
state, being -B(v) a unitary operator; finally, the rotation operator -R(p; v) 
that is function of the numerical vector p , has been introduced to represent 
the spin rotation produced by the operator g of eq.(2.4e). 

For completeness we also introduce the standard two components spinors w„ 
to represent the spin states. In this way the wave function corresponding 
I p, cr > is written as 

V'pa(q) =< q|p, o- >= wj(ci - p) (2.8) 
and the matrix elements of the boost operator of eq.(2.7) take the form 



< q, A = +1, //|S(v) |p, A = +1, a pP^^^]y^R^,{p; v)<5(q-p,(p; v)) 

' ^ (2.9) 
Here R^^{p;\) is the 2x2 matrix representation of the operator i?(p;v) 
acting in the space of the spinors w^.. 

We shall now generalize the procedure outlined above in order to include in 

the theory also the negative energy states. These states, as it is shown by 
the study of the Dirac equation and by the development of the field theories, 
are introduced for a consistent relativistic treatment of interacting particles. 
To this aim we have to replace, in eq.(2.4b), the positive energy Hamiltonian 
with 

/i = /i(A, p) = Ae(p) (2.10) 

where A represents here an operator with the eigenvalues A = +1 and A = 
— 1 for positive and negative energy, respectively. In consequence, one can 
inmediatly verify that the Poincare group commutation rules of eqs(2.3a- 
g) can be satisfied with the Hamiltonian of eq.(2.10) and by replacing, in 
eq.(2.4c) 

t = A(d + g) (2.11) 

By using the previous expression, the generalization of the finite boost for 
including the negative energy states is easily found: in eq.(2.5) the argument 
V must be replaced by Av, in agreement with the Lorentz transformation of 
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the three-momentum that was discussed before. For the generahzed boost 
operator we shall keep using the notation i?(v). 

In summary, for an Hamiltonian with negative eigenvalues (A = —1), also 
the boost generator must take a minus sign to give the correct commutation 
rules. The other generators, p and j of eq.(2.4a), remain unchanged. Having 
introduced for the free particle state the following ket |p. A, cr > we now 
need, to represent these states, a 4-component spinorial wave function of the 
following form 

<Aa(q) =< q|P, A, a >= u''iX)wJ{q - p) (2.12) 

that satisfies standard ortonormahty properties. We denote this representa- 
tion of the states as Poincare representation. In the previous equation we 
have introduced the following 4x2 components block spinors 

u^{+) = (l) (2.13a) 



«^(-) = (;) (2.136) 

for positive and negative energy, respectively. These spinors act, in eq.(2.12), 

on the standard two component spinor Wcr- In the Poincare representation 
the operator A is a 4 x 4 block matrix of the form 

One straightforwardly obtains the expressions for the Hamiltonian h and for 
the boost generator k replacing the previous expression of A in eqs.(2.10) and 
(2.11). Boosting the state |p. A, a > represents a generalization of eq.(2.7) 

S(v)|p, A, a >= [^M^^] V2i?(p; Av)|p,(p; Av), X,a> (2.15a) 

In consequence, the boosted wave function is of the form 



CA.(q) = [^]'^'u^mip; Av)^.5(q - p,(p; Av)) (2.156) 

where the equality q = p6(p; Av) given by the delta function has been used 
in the normalization factor. 
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For the matrix elements of the boost operator, eq.(2.9) is generahzed in the 
following way 



< q, A',/x|5(v)|p, A,(T >= 5x'x[ 



e(p&(p; Av ))ii 
^(P) 



^^^/xa(p; Av)5(q - P6(p; Av)) 



(2.15c) 

We now study the connection of the Poincarc representation with the one 
given by the solution of the standard Dirac equation. In particular this study 
is very useful in order to construct, in the following, invariant interaction 
operators with respect to boost transformations. We show that the Poincare 
representation is completely equivalent to the Dirac one. To this aim we recall 
that the solutions of the Dirac equation for a free particle, in the momentum 
space, with the same notations introduced before, have the following form 



V'g,,(q)=ii^(A,p)^,5(q-p) 



with 



^^''(+,P) 



u 



V,P) 



2e(p) 



2e(p) 



Ve(p) + 



m 



-\/e(pj+" 



e(p) + m 



(2.16) 



(2.17a) 



(2.176) 



In the previous equations the three Pauli matrices a have been introduced. 
Note that in the Dirac representation the spinors do depend on p (see the pre- 
vious equations), while in the Poincare one ( see eqs.(2.13a,b)) they do not. 
It is possible to pass from the former representation to the latter by means 
of the well-known, unitary but momentum dependent, Foldy-Wouthuysen 
(FW) transformation [25] given in the following equation 



e(p) + m 



nl/2 



+ 



e(pj — m 
26(p) 



nl/2 



>7) 



(2.18) 



26(p) J 

where we are introducing the Dirac matrices 7'* = (7°, 7) in the standard rep- 
resentation. With straightforward handling one verifies the following prop- 
erties of the FW transformation 



[/-i(p) = [/+(p) = C/(-p) 



(2.19) 
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The transformation for the spinors have the form 

[/(p)?i^(A, p)w„ = u^iX)tVa (2.20a) 

Note that, due to the unitarity of the FW transformation, the spinors u^{\, p) 
and u^{X) have the same normahzation to unity. Also, for the complete wave 
functions, defined in eq.(2.12), one has the relation 

C/(q)ii^(A, p)wj{q_ - p) = u''{X)wJ{q_ - p) (2.206) 

Applying the FW transformation as in the previous equation to the boosted 
Poincare wave function given in eq.(2.15b), the boost transformation of a 
Dirac wave function is easily found in the form 

Cp.(q) = ^^(q)<Ap.(q) = 

= f-^]'V(X, q)it:(p; Xv)wJ(q - p,(p; Av)) (2.21) 

where, analogously to eq.(2.15b), the equality q = P6(p;Av) given by the 
delta function has been used both in the normalization factor and in the 
argument of the Dirac spinor. 

On the other hand it is well known that in the Dirac theory the spinor boost 
is introduced in the form 

5^(v) = 55(v) = [1(.° + 1)]V2 + [1(.° - 1)]^/^^ - 

~l + ^(u7°7) (2.22) 

where the time component of the four- velocity, given in eq.(2.1d), has 
been used. Standard calculations show that 



Bn{v)u''{X, p)w, = [liP^ig^] V2^D(^^ p^(p. Av))it:(p; Av)^. (2.23) 

and, consequently, the boosted wave function of eq.(2.21) can be written as 
Cp.(q) = Bniy^^'^iX, p)wj{q - p,(p; Av)) (2.24) 
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For further developments, we note that eq.(2.23) can be simphfied by intro- 
ducing the covariantly normalized Dirac spinors 

«^^(A,p) = [^]V\^(A,p) (2.25) 
m 

satisfying the condition u^'-^{X, p)u^'^{X, p) = A. The Dirac boost for these 
spinors is 

5d(v)m^^(A, p)w, = M^^(A, (p6(p; Av)))i?(p; Av)w, (2.26) 

We recall that by using Dirac boosts and Dirac matrices one can construct 
Lorentz covariant operators. From standard algebra of the Dirac matrices 
one has 

B^(v)7°5z,(v) = 7° (2.27a) 

and 

Bz,(v)7V5d(v) = L^,(v)7°7'^ (2.276) 
In consequence, by using eq.(2.26) and eq. (2.27a), one obtains 

^+i?+(p; Av)u^^(A, (p,(p; Av)))«^^(A', (p,(p'; A'v)))i?(p'; AV)^., = 

= w+u'^^iX, p)ii^^(A', p')^;^' (2.28a) 
for the scalar matrix element. 

Also, by using eq.(2.26) and eq.(2.27b), one has 

w^R+ip; Xv)u^^{X, (p,(p; Av)))7''«^^(A', (p,(p'; AV)))i?(p'; XW)w,, = 

= L'^Mw^u^'^iX, p)7'^«^^(A', p')w,, (2.286) 

for the vector matrix element. Similar equations hold for the other Dirac co- 
variants, namely, the pseudoscalar, axial-vector and tensor matrix elements. 
Eqs.(2.28a,b), are very important for the construction of invariant interaction 
operators that will be done in the next section. 
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Prom eq.(2.26) one obtains the following useful expression for the spin rota- 
tion matrix 

R{p; v) = «^^(+, P6(p; v))i?(v)w^^'(+, p) (2.29) 

From the previous expression, with standard Dirac algebra one also finds the 
following relation 

i?+(p;v) = i?(p,(p;v);-v) (2.30) 
that will be used in the appendix to show the covariance of the interaction. 

3. Systems of relativistic interacting particles 

In this section we shall construct a Relativistic Hamiltonian Dynamical model 
for N interacting spin 1/2 particles also considering negative energy states. 
This objective will be achieved by defining the 10 total generators of the 
Poincarc group, denoted with the capital letters [G^], in terms of single par- 
ticle operators. Obviously, the total generators must satisfy the same com- 
mutation rules given in eqs.(2.3a-g) for the single particle generators. If the 
interaction were not present, one could easily define the total generators as 
the sum of the single particle ones: 

TV 

G' = T.9/ (3.1) 

1=1 

automatically satisfying the commutation rules. 

As anticipated in the introduction, the way in which the interaction is intro- 
duced makes the difference among various models of RHD. For this problem, 
that has been mainly faced considering positive energy states, different solu- 
tions have been proposed as explained in subsect.1.1. 

In the present work we want to keep using standard Lorentz transformations 
of the four-vectors (canonical boosts), as given in eqs.(2.1 a-c) and (2.2a,b), 
so we can only choose IF RHD (i) or PF RHD (ii). 

(i) As for the IF RHD, the interaction is added to the sum of the free Hamil- 
tonians but not to the momenta. In consequence, an interaction operator 
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must be also added to the sum of the free boost operators in order to satify 
the Poincare algebra commutation rules. This method, that is also adopted 
for the quantization of the relativistic field theories, has been widely used 
to introduce relativistic corrections [2] both to the Hamiltonians of bound 
systems and to the operators that describe the interaction of these systems 
with external electromagnetic fields, significantly improving, for the quark 
models, the reproduction of the experimental data [26], specially for the low 
energy observables. 

The difficulty of this approach consists in finding the exact expression of the 
interaction dependent operator that modifies the boost generator. For this 
reason we do not follow this method in the present work. 

(ii) In the PF RHD, the interaction modifies both the Hamiltonian and the 
total momentum of the system, leaving the boost free of the interaction 
[5-7]. As shown in the following, the form of the interacting four-momentum 
operator can be directly determined. 

A relevant consequence of the properties of PF RHD is that it is possible to 
study the dynamics of the composite system in terms of explicitly covariant 
integro-differencial wave equations that will be derived in the next section. 

We revise the Point Form procedure considering the possibility of introducing 

also negative energy states. 

Given a (bound) system of N interacting constituent particles, it is possible 
to observe this system both in its rest reference frame (RF) and in a generic 
reference frame (GF). It is convenient to introduce the observable quantity 

= V) 

that represents the four- velocity of the RF measured from a GF. It means 
that V is the parameter that, inserted in eqs.(2.2a,b), allows to transform the 
momenta observed in the RF of the system to the corresponding quantities 
in the GF. The relation between and V is the same as in eq.(2.1d). 
Furthermore, for a system of mass M, considering as a, classical quantity, 
one has 

E P 

={—,—) (3.2) 



where E = \/P^ + M'^ and P respectively represent its energy and three- 
momentum measured in a GF. 
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The procedure to construct the generators of the Poincare group requires 
to define V as an operator, that IS clS db dynamical variable of the system. 
This definition will be given in eq.(3.7). As first step we introduce the RF 
four-momentum of the i-th particle 

pr(A.) = (A,6(p*),p*) (3.3) 

where the asterisk denotes the quantities observed in the RF. The sum of 
these four-momenta over the N constituents, by definition of the RF (that is 
also called zero momentum frame), is given by the following equation 

JV N 

EPr(A.) = (EMP*) = M^,0) (3.4) 

1=1 i=i 

where we have also introduced Mp that represents the free mass operator of 
the system. By applying the Lorentz transformation of eq.(2.2b) as function 
of the parameter V to the p*^{\i), also using eq.(3.4), one can write the sum 
of the four-momenta of the particles in a GF as 

N 

EP^iK) = V^Mp (3.5) 

i=l 

If Mp is nonvanishing one can solve the previous equation with respect to 
V^; then by writing Mp in terms of the Pi{Xi), one can express as a 
function the (Aj), or, more precisely, of the and Aj, that are chosen as 
dynamical variables of the relativistic model. 

The condition Mp 7^ is fulfilled by the states in which all the particles 
have the same energy sign, that is Aj = A for i = 1, A^. These states will 
be denoted as dynamical states. On the other hand, the states in which the 
particles do not have all the same energy sign can give a vanishing value of 
Mp, not allowing for a definition of in terms of the particle momenta. 
For this reason, these states will be treated separately as auxiliary states. 
For the dynamical states the expression of Mp as a function of the momenta 
in a GF is the following 

1 

(3.6) 

where we have introduced the collective shorthand notation {p} = pi, pjv. 
Analogously we introduce the notation {A} = Ai, Ajv. For the dynamical 
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N 



Mp^Mp{A,{p})=A 



states one has {A} = Ai = A2 = ... = A^r = A. In consequence, we can also 
write 

N 

V^{{X}^A,{p}) = [M^(A,{p})]-^EPr(A) (3-7) 

i=l 

Let us note that the observable four- vector V^, as given in the previous 
expression, for both A = +1 and A = — 1, transforms in the same way as a 
positive energy four-momentum, that is replacing p with V, e(p) with y^(V) 
and setting A = +1 in eq.(2.2b). In this way we introduce 

F,° = \/,°(V;v) = l-°(V,(V;v)) (3.8a) 

V6 = V6(V;v) (3.86) 

This result, that is consistent with the definitions of eqs.(3.2a-c), can be 
easily derived by transforming, with the help of eq.(2.2b), the p^(A) that 
appear in eq.(3.7). 

We shall now choose the complete set of commuting operators that will be 
used for the quantum-mechanical description of the system . To this aim we 
note that, due to its definition in eq.(3.7), the operator commutes with 
the momenta of all the particles. In consequence, it is possible to choose the 
following operators: the three-momenta of N — 1 particles pi, Pn-i = {q}, 
the spatial components of the four-velocity V, the energy signs {A}, and, 
finally, the spin projections on the z axis ai, ...,aN = {c}. 
With this set of commuting operators, the representation states that will be 
used to write down the wave functions of the model, are of the form 

|V'.>=|{q},V,{A},M> (3.9) 

Their normalization is 

< ^I^M >=< {q},V,{A},M|{q'},V',{A'},{<7'} >= 

= S'ip^ - p[)...6%pN-i - p'N-i)S'{y - V')5wK}5{A}{A'} (3.10) 

We choose these representation states in order to derive in a simple way the 
manifestly covariant wave equation of the model and, in turn, to obtain a 
dynamically conserved current, as it will be studied in subsequent works. 
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A different type of representation states, currently denoted as velocity states 
can be advantageously used to study the relativistic bound state wave func- 
tions. In the velocity states the spatial variables are represented by V and 
by the N (not indepedendent) rest frame momenta {p*} or better by the 
N-1 (independent) Jacobi momenta {k}. As shown in ref. [5], the Lorentz 
transformation of these states is given by the standard boost of V, as in 
eq.(3.8b), and by a Wigner rotation for the {p*} or for the {k}. 

In general, we point out that, in our relativistic model, as it will be shown 
in the following, the total three-momentum P is interaction dependent, so it 
does not commute with the three-momenta of the constituent particles pj. 
For this reason P cannot be diagonalized simultaneously with them and, as 
discussed before, V is conveniently chosen from the beginning. 

We note that the momentum of the N*^ particle, when it appears in the 
calculations, can be expressed as function of {q}, V and {A}. To this aim, 
we firstly introduce the four-vector 

N-1 

Q'=T.Pfi>^^) (3-11) 

i=l 

We also recall the standard relation 

A.e(p*) = V,p^(X,) (3.12) 

Then, we write eq.(3.5) with the definition of Mp given in eq.(3.4) in the 
form 

N 

Q" = -p'nM + V''Y.>^,e{p*) (3.13a) 

i=l 

Squaring both sides, with the help of eq. (3.12), we find 

e{p^*) = e;,({A}, {q}, V) = [{Q^V,)' + m'- Q^Q, ' 
Then 

({A}, {q}, V) = V^Q^ + AAre^({A}, {q}, V) (3.13c) 
and finally, by means of eq.(3.5) 

P%m, {q}, V) = -Q^ + V^Mj.{{X}, {q}, V) (3.13ci) 



(3.136) 
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The previous expression can be also used for the auxihary states. In fact, 
when Mf{{X}, {q}, V) = 0, one has 



P^^({A},{q},V) = -g'' 

In the case of dynamical states, for the developments of the next section it 
is convenient to introduce 

Mp(A = -1, {q}, V) = -Mf{A = +1, {q}, V) = -Mf{{ci}, V) (3.13e) 

Explicit expressions for the positive free mass Mi?({q}, V) will be given in 

eqs.(4.8) and (4.16a,b) for the two and three-body case, respectively. 

We can now take advantage of eq.(3.2) to define the total momentum operator 

as 

pi^^MV^ (3.14a) 
Here M represents the invariant mass operator of the model, defined as 

M ^ MF{{X},{q},y) + W (3.146) 

where W represents the Lorentz invariant interaction operator. Our pro- 
cedure for introducing the interaction represents the generalization of the 
Bakamjian-Thomas construction [5,27] to a theory with negative energy states. 
In order to obtain a Lorentz invariant operator, we require the following 
commutation rule of W with the boost generator 

[W, T] = (3.15) 

We also require 

[W,V] = (3.16) 

From the previous equation, recalling that = [V^ -|- 1]^^^ one straightfor- 
wardly has 

[W, y°] = (3.17) 

ensuring that no commutation problem arises when defining the total mo- 
mentum P*^ in eq.(3.14a). Furthermore, from the definition of eq.(3.14a) and 
the requirement of eqs.(3.16) and (3.17) one has 

[P'^,M]=0 (3.18) 

ensuring the space and time translational invariance of the model. 
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In order to satify eqs.(3.15) and (3.16) we choose the interaction operator in 
the form 

^= E ^{A}{A'} (3.19a) 

{A}{A'} 

with 

W{XHy}- E / ^^{qK{qVV F({q}, {q'}) 

{a}{a'} 

< {ci},Y,{X},{a}\W^\W,y,{X'},{a'} > 

|{q}, V, {A}, M X {q'}, V, {A'}, {a'}\ (3.1%) 

where the integration rf^{q} simbohzes the integration over the momenta 
Pi,...., Pat-i and analogously for the primed variables; F({q},{q'}) is a 
function of the momenta, that will be explicitly given in eq.(3.20), that en- 
sures the Lorentz covariance of the total operator. Also, we introduce the 
manifestly covariant interaction amplitude , in the form 

< {q}, V, {A}, {a}\W^\W, V, {A'}, {a'} >= 

N 

E E Kf({q},{q'},{A},{A'},v) 

K i>j=l 

w+{.yu^^{{\h {pmfym, {p'})^{.o (3-19C) 
in the previous equation 

^#({q},{q'},{A},{A'},v) 

is a Lorentz invariant function, i.e. depending on the scalar products of the 
four-momenta, that expresses the spatial part of the two-body interaction of 
the model; 
rfj can represent: 

a) the product of the covariant Dirac matrices (introduced in the previous 
section) for the particle i and j; more explicitly, for = 1 one has a scalar 
interaction with lilj , being 7j the identity Dirac matrix of the i*'* particle; 
for K = 2 one has a vector interaction with '^i'yj'gfj.u and so on for the 
pseudoscalar, axial-vector and tensor interactions; but also, for X > 5 
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b) terms containing Lorentz invariant products of the Dirac matrices with 
the four-momenta of the — 1 particles, or with V^, that is V]u7f , as we 
shall see in the next section for a specific model. 

Furthermore, the notation u^'-^{{\}, {p})w^uj represents the direct product 
of the Dirac spinors (see eq.(2.25)), for all the particles. Here and in the 
following, p'^ and the corresponding primed quantity are given by eq.(3.13d). 

The structure of eq.(3.19b) immediately shows that the interaction opera- 
tor W satisfies the commutation rule of eq.(3.16). In fact, the interaction 
amplitude of eq.(3.19c) depends on V but the operator W of eq.(3.19b) has 
vanishing matrix elements between states with different values of V. Fur- 
thermore, the Lorentz invariance of the interaction, that is the commutation 
of the interaction operator with the boost generator, expressed by eq.(3.15), 
is ensured by the invariance of the Dirac spinor matrix elements that appear 
in eq.(3.19c) taking, in eq.(3.19b), the function F({q}, {q'}) in the form 

i^({q},{q'}) = [<{q})e({q'})]-^/' (3.20) 

The details of the demonstration are given in the appendix. According to 
eq.(A.2), here and in the following subsection ( in particular when considering 
the covariant wave equations) the three-momentum of each particle must be 
transformed according to its energy sign as in eq.(2.2b). 

We have now all the elements to prove that the generators [G^] of our rela- 
tivistic model satisfy the Poincare algebra commutation rules of eqs.(2.3a-f). 
The only generators that contain the interaction are = H and P de- 
fined by means of eqs.(3.14a,b). The other ones, being free of the interac- 
tion, are given by the sum of the single particle generators. For this reason 
eqs.(2.3b),(2.3d) and (2.3f) are automatically satisfied. 
From the definition of {P^ = H,P) given in eqs.(3.14a,b) one immediately 
obtains the first two relations of eq.(2.3a). 

From the rotational invariance of V^, Mp and W, one obtains the last relation 
of eq.(2.3a). 

Taking also into account the vector character of V, one verifies eq.(2.3c). 

Finally, considering the definition of V in terms of single particle operators 
as given in eq.(3.7), with [T, Mp] = and eq.(A.l) for the interaction term, 
one simply derives eqs.(2.3e,g) completing the verification of the Poincare 
invariance of the model. 
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4. The wave equation of the model 

In this section we study explicitly the eigenvalue wave equation of the model 
according to the properties discussed in the previous section. In the sub- 
sect. 4.1 we shall analyze the general structure of the equation, higlighting 
its manifest covariance, while in the subsect.4.2, by using some techniques 
developed by RWE, we shall introduce the auxiliary states and the quadratic 
terms of the quasipotential for two and three-body systems. 

4.1 General structure of the wave equation 

The mass eigenvalue equation, by means of the mass operator definition of 
eq.(3.14b), may be written in the general form 

D(M,{q},A,V)|* >= > (4.1a) 

with 

D(M, {q}. A, V) = M - Mp({q}, A, V) (4.16) 

Here we use the free mass operator Mp introduced in eq.(3.13c) only for 
the dynamical states; the structure of the interaction operator W has been 
given in eqs.(3.19a-c) and (3.20); the operator D(M, {q},A, V) has been 
introduced here only to simplify the comparison with RWE models. 

Considering only dynamical states, one has 

I* |^,A = +1,V > +|*,A = -1,V > (4.2) 

Projecting eq.(4.1) onto the states defined in eq.(3.9) gives the following set 
of coupled equations 

[M - AM^({q}, V)]*({q}, V, A, {a}) - 
Y: [ d'W < {q}, V,A,M|iy|{q'}, V,A',M > *({q'}, V, A', {a'}) 

A' {a'} 

(4.3) 

where also the definition of the positive free mass operator given in eq.(3.13e) 
has been used. 
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In order to find a numerical solution, the previous equation can be con- 
veniently written in the RF of the bound system, then the obtained wave 
function is standardly boosted to any GF. 

In the RF the sum of the N three-momenta of the particles gives zero, as 
shown in eq.(3.4). This allows to introduce the N — 1 independent Jacobi 
momenta {k} that are used to study in a clear way the symmetries of the 
spatial part of the wave function. In this way cq.(4.3) takes the following 
form (for brevity, here and in the remainder of the paper, we do not write 
the RF eigenvalue V = 0): 

[M-AMH{k})]*({k},A,M) = 

E / d'{k'} < {k},A,{a}\Wm,A',{a'} > ^({k'}. A', {a'}) (4.4) 

A'{a'} 

We now turn to write eq.(4.3) in a manifestly covariant way. By introducing 
$({q}, V, A, {a}) = [e({q})]^/^vl;({q}, V, A, {a}) (4.5) 

and with the definition of the manifestly covariant interaction amplitude of 
eq.(3.19c), eq.(4.3) is written as 

[M - AMp({q}, V)]$({q}, V, A, {a}) = 
^ /^^<{q},V,A,M|VF^|{q'},V,A',K}> 

A'{a'} ^\I^J) 

$({q'},V,A',K}) (4.6) 

Note that the integrations over the — 1 particle momenta is performed in 
a covariant way by means of the factor e({q'}) in the denominator of the 
r.h.s. of the previous equation. A similar structure of covariant integration 
was firstly used in the RWE proposed by Gross [13]. 

We point out that, as manifest covariance explicitly shows, PF RHD allows 
to boost in an unambigous way the wave function of the model for calculating 
physical observables. 

The covariant interaction amplitude of eq.(4.6) should be determined in or- 
der to reproduce, with the best possible approximation, the dynamics of 
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the underlying field theory. In particular, the inclusion of the z-graphs will 
be performed in the next subsect.4.2 obtaining the amplitudes of W^j^j^ in 
eqs. (4.13b) and (4.17b) for the two and three-body case, respectively. 

For completeness, note that the standard form of PF RHD could be recov- 
ered by completely excluding in eq.(4.6) the negative energy states, that is 
taking only the first term in eq.(4.2). Obviously, this choice does not violate 
the relativistic invariance of the model. In this way, as anticipated in the 
introduction, the so called Blackenbecler-Sugar wave equation [9] is obtained 
without difficulties and the Poincare algebra commutation rules are satisfied 
by generators that only act onto positive energy states. 

Another possible choice is to include the states with A = — 1 and to use an 
interaction operator without z-graphs. It corresponds to the Breit equation. 

In order to make a comparison with other relativistic models, we consider 
explicitly the case of a two-body system, that is particularly relevant in nuclear 

physics for the study of the deuteron. 

In this case the only Jacobi variable is the relative momentum k = = — P2 
and eq.(4.4) for the RF takes the form 

[M-2Ae(k)]*(k,A,M) = 
^ / d^k' < k,A,{a}\W\k',A',{a'} > *(k'. A', {a'}) (4.7a) 

A'{a'} 

The manifestly covariant form of the previous equation (to be used in a GF), 
derived directly from eq.(4.6), is 

[M - AMf(pi, V)]$(pi, V, A, {a}) = 

E / ^ < Pi. V, A, Mliy^lp/, V, A', {a'} > $(p/, V, A', {a'}) 

(4.76) 

where, by means of eqs.(3.13a-e), we can express the positive free mass op- 
erator as 

Mp{p^,V)^2p^,V^^2e{k) (4.8) 

In eqs.(4.7a,b) the two particles have the same energy sign as in the Breit 
equation that is obtained by reducing to a three-dimensional form the Bethe- 
Salpeter equation in the case of a static interaction between the particles. In 
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particular, in ref. [14], the Coulomb interaction was considered to derive the 
Breit wave equation in the RF of the bound system. 

To reproduce that result in our model, the following covariant interaction 
amplitude must be inserted in eq.(4.7b) 

< Pi, V, {A}, {a}\WEM, V, {A'}, {a'} >= 

CBr(Al, A2, X[, X2)UBr{y, Pi, P'l) 

w+,,«^^(Ai, pi)w+^^u^^{X2, P2) 

( W)(K72)«''''(A'i, p;)«;.;«^^(A'2, p'^w,'^ (4.9a) 

with 

CBr({A},{A'})=Ai5A,A,5AiA^ (4-96) 

and the covariant function 



C/b.(V,Pi,p;,{A},{A'}) = -m^^[|\^X(AO Kpi' "(KW^' 

(4.9c) 

In the RF the two-body Breit equation is obtained in the standard form of 
eq.(4.7a) with 

<k,A,{a}\WBr\k',A',{a'} >= 
-m'^,CBr{{X},m) 
<n^+(Ai, k)w+y+{X2, -k) _^j^,^, h^(A;, k')w,'y{X'„ -kO^.^ (4.10) 

Summarizing, we note that our procedure based on the commutation rules 
of the Poincare algebra, naturally introduces the covariant general expres- 
sions of eqs.(4.6); eqs.(4.7a,b) arc simply obtained by specializing eq.(4.6) 
for a two-body system. On the other hand, the form of the interaction is 
not determined by Poincare algebra and must be chosen according to a spe- 
cific dynamical model for the bound system. As an example, in the case 
of the Breit equation discussed above, the instantaneous approximation is 
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performed to sum up in a three-dimensional form the ladder, or uncrossed, 
Feynman graphs introduced in the Bethe-Salpeter equation. 

As for the states with different energy signs, we can analyze in more detail the 
point shown in the previous section. For these states, the free mass operator 
Mi?({q}, {A}, V) can be, in general, vanishing and in particular, for two- 
body systems with {A} = (+, — ) and {A} = (— , +) it is always vanishing, 
being e(p^) = e(— Pg) = e(k). In consequence, it is not possible to define, by 
means of eq.(3.7), the dynamical variable V in terms of the three- momenta 
of the particles . It means that Poincare invariance, in the form discussed 
in sect. 3, only allows to introduce the dynamical states of eq.(4.2). 
Such prohibition given by the Poincare algebra naturally leads to the dy- 
namical consequence of excluding an unphysical pole at M = in the Green 
function of the model, avoiding the continuum dissolution or cockroach nest 
disease [11,12]. 

4.2 Auxiliary states and quadratic terms 

We first consider the two-body case (i), then we generalize the model to the 
three-body case (ii). 

(i) Two-body case. 

The difficulty of the M = pole was overcome in the Mandelzweig- Wallace 
(MW) model, by performing, for two particle bound systems, a three-dimen- 
sional reduction of crossed and uncrossed Feynman photon exchange graphs. 
This procedure, also denoted as equal time reduction, making use of the 
eikonal approximation for the z-graphs, is more accurate than the one of the 
Breit equation and correctly reproduces the scattering T matrix up to the 
terms of second order in the interaction operator [11,15,17]. 
One has a wave equation written by means of the inverse of the Green func- 
tion that in the RF (V = 0) takes the form 

DMwiM, k, {A}) = 6x,x,M + rf(Ai, A2)2e(k) (4.11) 

with d{+l, +1) = -1 and d{+l, -1) = d{-l, +1) = d{-l, -1) = +1. 

The MW equation can be used with different kinds of interaction operators. 

Without entering here into details, one can use an operator of the same type 



32 



as that given in eqs.(4.9a-c), replacing CbtH^}, {^'}) ^^^^ 

Cuwm, {A'}) = 1 - 2dx,x,Sx',x',Sx,-,5y^-i (4.12) 

allowing for matrix elements with the states that have different energy signs. 
However, in that model, it is not possible to write the wave equation by 
defining the mass operator, as required by the Poincare group commutation 
rules in the framework of the Relativistic Hamiltonian Dynamics. 

Considering as starting point a two-body system, we propose to keep the def- 
inition of the mass operator given in the present work and to use eqs.(4.7a,b) 
as the dynamical wave equation of the model. We include the states with 
different energy signs, that represent the z- graphs, by means of a re-definition 
of the interaction operator and by the explicit introduction of the auxiliary 
states. 

We introduce the following effective interaction operator that contains the 
z-graphs in the second quadratic term, in the form 

<k,A,{a}\W,ff\k',A',{a'} >= 

<k,A,{a}\WMw\\^',A\{a}> + 
Y: [ d'k" < k. A, {a}\WMw\i^", {A"}, {a"} > 

{A"}{a"} 

< k", {A"}, {a"}\WMw\i^', A', W} > (4.13a) 

where, in the last term, the states with mixed energy signs are represented 
by {A"} = (-1-1, —1), (—1,-1-1). The corresponding manifestly covariant am- 
plitude is: 

< pi, V, A, {a}\W^ff\p,', V, A', {a'} >= 

< Pi, V, A, {a}\WZ^\p,', V, A', {(j'} > + 

^ /^<Pi,V,A,M|li^l^^^|p/',V,{A"},{a"}> 

[Mf(Pi", V)]-^ < p^",y,{X"},{a"}\W^^\p,',V,A',{a'} > (4.136) 
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where we have used the factor of eq.(3.20) and we have expressed the de- 
nominator of eq.(4.13a) by means of cq.(4.8). 

For the calculation of the bound state properties, the effective interaction 
matrix elements of eqs.(4.13a) and (4.13b) must be inserted in the eigenvalue 
equations eqs.(4.7a) and (4.7b), respectively. 

Also, assuming that the states with equal energy signs represent the dominant 
contributions in the MW equation, we can define, in our model, the auxiliary 
states as 

*(k,{A},M) = 

E /d3k'<k,{A},M|VFM^|k',A',{a'}>*(k',A',{a'}) (4.14a) 

The corresponding covariant definition is 

$(pi,V,{A},M) = [M^(pi,V)]-^ 

^ /^4<Pi,V,{A},M|VF^^|p;,V,A',{a'}>$(p;,V,A',{a'}) 
A'{.'}-^ e(pi) 

(4.14&) 

Considering a perturbative expansion up to the second order in the inter- 
action, one has the same result that is given by the MW equation. In this 
sense, our model is dynamically equivalent to that by MW. 
Furthermore, as explained in subsect.1.2, the explicit introduction of the 
auxiliary states, even if they do not participate directly in the dynamics, can 
be very important for the calculation of the electroweak interaction matrix 
elements of the bound system. 

{a) Three-body case. 

We now propose the structure of a possible generalization of the model for 
a bound system with = 3, that is relevant for baryonic quark models, 
including the states with mixed energy signs. The same procedure can be 
straightforwardly generahzed to the case of A" > 3 particles. 
First, we specialize eqs.(4.4) and (4.6) to the case of a three-body system, 
obtaining in the RF 

[M - AM^(p„ Pa)]^'(Pp, Pa, A, {a}) = 
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Y: [ d%d'p', < p„ p„ A, {a}\W\p'^, p'„ A', {a'} > *(p;, p',, A', {a'}) 

A'{a'} 

(4.15a) 

where we have introduced the standard RF three-body Jacobi momenta 
Pp = - P2) and Pa = ^Jl{pl + vl). 

The manifestly covariant form of the previous equation (to be used in a GF) , 
is derived directly from eq.(4.6) in the form 



E 



[M - AM^(pi, p2, V)]<l>(pi, P2, V, A, {a}) = 
d^p[ ^'Pi < p,, p„ V, A, {a}|l^^|p;, p^, V, A', {a'} > 



^p[,p'„V,A',{a'}) 
The positive free mass operator in eq.(4.15a) is 

Mf{Pp, Pa) = e(pD + e{pl) + e{pl) = 



(4.156) 



:PA + y-Pp| +m 



+ 



1/2 



+ 



gPA-y-Pp 



1/2 



m 



1/2 



(4.16a) 



With the help of the cqs.(3.11), (3. 12), (3. 13b) and (3.13e), the positive free 
mass operator in eq.(4.15b) is 



Mf(Pi, P2, V) = (p^+p^) V;+ [((K + p'2)V,y + m'- {p1 + pi^) {p\ + p^,)g,, 

(4.166) 

For the effective interaction, in order to take into account the effects of the 
states with mixed energy signs, we use the following expression 



1/2 



< Pp, Pa, a, {a}\Weff\p'p, Pa, A', {a'} >= 
< Pp, Pa, a, {a}\W\p'p, p^. A', {a'} > + 

Y: [ d'p^pl < p„ Pa, a, {a}\W\p';, pi {A"}, {a"} > 



{X"}W'} 
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. I rxn\ < Pp,Pa,{A"},K}I^Ip>'a,A',{^'} > (4.17a) 

For the corresponding covariant amplitude we have 

< pi, p2, V, A, MlW^g^lp;, p'2, V, A', {a'} >= 

< pi, p2, V, A, {a}\W^\p[, p'2, V, A', {a'} > + 

^ / I'^M I'^M < P^' ^' ^' {^ll^^lPi, P2, V, {A-}, {a-} > 

g(p;.p;!v.{A»}) < {""il^^lP'.. P'- V. ^' > 

(4.176) 

As in the two-body case the second quadratic term represents the z-graphs 
of the model. 

The auxiliary states are expressed in general form as 

*(Pp,PA,{A},M) ^ 



i^(Pp,PA,{A}) 

Y: f d%d'p\ < p„ p,, {A}, {a}\W\p'^, p'„ A', {a'} > *(p;, p'„ A', {a'}) 

A'{(t'} 

(4.18a) 

the covariant form is 

$(pi,P2,V,{A},M) = 



5(pi,P2,V,{A}) 
E / <Pi,P2,{A},Mr^|pi,p'2,V,A',K}> 

$(p;,p'2,V,A',{<7'}) (4.186) 

where the notation {A} refers, as before, to the states with mixed energy 
signs. The specific form of the interaction amplitude should be fixed on the 
base of dynamical considerations that go beyond the scope of the present 
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work that essentially concerns Poincare invariance. For the invariant, non- 
vanishing quantity S(pp,pA,{A}) or i?(pi, p2, V, {A}) in a previous study 
devoted to three body systems [28] , we proposed a model in which they take 
the form 



S(p„pA,+l,+l,-l) = 2 



' 2 I 2 



1/2 



(4.19a) 



S(pi, p2, V, +1, +1, -1) = 2 \Mf{pi, P2, V) - VM + P2) 



(4.196) 



and 



B{pp,p^, +1,-1,-1) = 




m 



1/2 



Mf{Pp,Px) 

(4.20a) 



B{p^, p2, V, +1, -1, -1) = - ^i^(Pi> P2, V) (4.20&) 

The results for other values of {A} are determined by symmetry reasons. 
This model consists in a generahzation of the MW equation to three body 
systems that gives the correct nonrelativistic limit. Further investigations 
must be developed on the dynamics of the model in order to determine the 
relevant Feynman graphs (and their three-dimensional reduction) for a three- 
body system. 

We conclude remarking that Poincare invariance, implemented by means of 
the Relativistic Hamiltonian Dynamics, offers a reliable framework to study 
in a covariant way the strongly interacting bound systems, also including 
negative energy states. The wave functions of the model can be exactly 
boosted to any reference frame. Manifestly covariant wave equations can be 
written avoiding some pathologies encountered in other models. 
Numerical investigations are required to study the effects of negative energy 
states on the electroweak form factors. 

At a more fundamental level, it would be of great interest to study in detail 
the connection between the Relativistic Hamiltonian Dynamics and the for- 
malism of an underlying field theory, that, for the quark models, is assumed 
to be the QCD. 
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Appendix 

Instead of proving directly eq.(3.15), we study the finite boost transformation 

B(v)WB+(v) = W (Al) 

that, expanded at the first order in v, is equivalent to cq.(3.15). Furthermore, 
in this way we also learn how the quantities that appear in the model are 
transformed under finite boost. The procedure can be summarized in the 
following steps. 

(i) We write down in the next equation the explicit expression of the action of 
the boost operator onto the representation states of eq.(3.9). By considering 
the boost of eq.(2.15a) for a single particle, one obtains for the total boost 

5(v)|{q},V,{A},W 

i?({p},{A};v)|{q,}, V6,{A},M > (A2) 

In the previous expression , {qh} = {q6}({q}, {A}; v) denotes the set of the 
standardly boosted momenta for the particles 1, — 1 as functions of the 

unboosted momenta, of the energy signs and of the parameter of the Lorentz 
transformation, respectively; e({q5}) represents the product of the corre- 
sponding boosted absolute values of the energies; recalling that V must be 
considered as a function of the particle momenta as shown in eq.(3.7), when 
a boost operator is applied to eigenstate of V, one obtains an eigenstate of 
Vb = Yb{Y] v). Its expression and that of have been given in eqs.(3.8a,b); 
finally, -R({p},{A};v) represents the product of the spin rotation matrices 
for all the particles. 

(n) We calculate B{v)WB^{v) by using the operatorial relation of eq.(A.2) 

onto the interaction defined in cqs.(3.19a,b) and then replace the integra- 
tion variables {q} with {q;,} and analogously for the primed momenta; we 
also replace V with V^. These change of variables introduce the following 
transformation factors 

d'id} = 1^ rf^{q.} (A.3) 



<{q5})n° 
e({q})^° 
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and the analogous one for the primed momenta; also 

d'y = ^ d'Y, (AA) 

At this point wc rc-define the integration variables eliminating the index b 
in {qb}, {qb'} and V(,, but now the arguments {q}, {q'} and V appearing 
in eq.(3.19b) must be expressed in terms of inverse Lorentz transformations 
(obtained with the boost parameter — v) of the new integration variables. 
The result is 

^(v)W^{A}{A'}i?+(v) = E [ d'{ci}d'{c^}d'y 

WhW'} 



E E V;f({q},{q'},{A},{A'},V) 

K i>j=l 



^({P."},{A};-v)|{q},V,{A},M> 
<W},y,{X'},{a'}\R-'{{p'f},{X'};-v) 



e{{ck})e{H-} ) 

e({q})<{q'}) 



nl/2 



iA.5) 



where the shorthand notation {p^ } = {P6}({P6}, {X}; — v) and the analogous 
for the other momentum variables have been introduced. 



(iii) By considering the sum over the complete set of the spin variables, the 
rotation operators can act on the spinors and w^^j . Then one uses 

for these operators the property of eq.(2.30) with the boost parameter equal 
to —V and the transformation of the Dirac covariant matrix elements of 
eqs.(2.28a) and (2.28b). 

Finally, by using the function defined in eq.(3.20) one verifies eq.(A.l) for 
each W{x}{x'} 
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